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SUMMARY 


I  ' 

For  a  nurr.ber  of  turbulent  flov;  problems,  a  first-order  closure 
or  eddy  transport  approach  to  the  modellnp;  of  the  turbulence  is 
not  sufficiently  accurate.  It  is  necessary  that  a  more  powerful 
apcrcach  be  used  to  study  the  dynamics  of  the  turbulence.  Flows 
involving  chemical  reactions,  such  as  those  in  a  chemical  laser, 
are  good  examples  of  flows  requiring  a  more  powerful  approach.  The 
nixing  in  HF  and  DF  chemical  lasers  takes  place  under  conditions  of 
very  large  heat  release.  The  coupling  between  the  heat  release  and 
the  turbulence  is  an  Important  feature  of  the  flow.  A.R.A.P.  has 
been  fund^^tm-d-er- APO-SR  Contracts  F4 4620-73-U-0027  and  P44R20-73-C- 
092^''^to  employ  this  basic  method  to  investigate  problems  connected 
with  chemical  lasers  of  particular  interest  to  the  Air  Force  V.'eapcns 
Laboratory.  These  studies  have  resulted  in  considerable  insight 
into  the  complex  processes  occurring  in  laser  cavities.'^  These 
studies  have  been  documented  in  References  1,  2,  and  3- 

^•This  report  describes  recent  results  of  studies  using 

the  reacting  shear  layer  (RSL)  computer  program  and  discusses  the 
proposed  procedure  for  handling  multi-step  chemical  reactions  in  a 
second-order  closure  computation  of  turbulent  reacting  flows. ^ 
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I  .  TMTRODUCTIOIJ 


A  second-order  closure  prop;ram  for  the  investir^at ion  of  tur¬ 
bulent ,  compressible,  chemically  reactinpr  flowfields  is  under 
development  at  A.R.A.P.  Second-order  closure  metliods  have  provided 
an  understand  j  nr  of  the  physics  of  turbulent  flov/s  which  could  not 
be  provided  by  classical  first-order  closure  techniques  such  as 
the  eddy  viscosity  method.  We  have  successfully  applied  these 
techniques  to  a  variety  of  turbulent  flov;  problems  includinp;  incom- 
pressible  and  compressible  shear  layers/'  the  behavior  of  the 
planetary  boundary  layers,  and  the  distiersal  of  chemically  reacting 
pollutants  in  a  turbulent  atmosphere.' 

A.R.A.P,  is  funded  under  AFOSR  Contract  ho.  ^620-7 j-C-0026 

to  develop  an  invariant  second-order  closure  model  capable  of  investi¬ 
gating  the  mixing  and  reactions  in  HP  and  DF  chem.lcal  lasers.  The 
mixing  in  these  devices  takes  place  under  conditions  of  large  local 
heat  release  due  to  the  highly  exothermic  chemical  reaction.  The 
interaction  between  the  heat  release  and  the  turbulence  is  an  impor¬ 
tant  feature  of  the  flow. 

The  set  of  second-order  closure  equations  governing  the  mixing 
and  reaction  of  compressible,  turbulent,  multi-component  reacting 
shear  layers  has  been  modeled  and  programmed  for  solution  in  our 
R.SL  program..  This  interim  scientific  report  describes  some  recent 
results  of  studies  using  this  program  and  discusses  the  proposed 
procedure  for  extending  the  program  for  studying  multi-step  chemical 
reactions  in  a  second-order  closure  computation. 
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II,  REACTIMli  5’>I!I;AR  LAYER  (RRL)  F^^onRAM 


The  reacting  shear  layer  program  has  been  completely  cFiecked 

out  ever  the  past  six  months.  Program  predictions  for  nonreacting 

6  T 

flows  ’  (two-dimensional  wake,  heated  wake,  heated  planar  jet) 
are  in  good  agreement  with  experimental  measurements.  Comparison 
of  program  predictions  with  experimental  data  for  reacting  flov/s 
is  now  underway.  The  RSL  program  has  been  used  for  computations 
of  the  HP  chemical  laser  flov;field,  hydrogen-air  premixed  combustion 
and  diffusion  flames  and  propane-air  premixed  combustion  and  dif¬ 
fusion  flames.  The  calculations  have  been  performed  using  two  dif¬ 
ferent  models  for  the  scalar  correlations.  One  model  (used  during 
the  model  development  stage)  was  to  set  all  third-order  (and  higher- 
order)  scalar  correlations  to  sero.  The  other  model  is  the  A.R.A.P. 
’’typical  eddy”  box  model.  Preliminary  results  for  exothermic  react- 
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ing  flows  suggest  that  there  is  only  a  small  difference  between  the 
calculations  using  these  two  models.  More  extensive  tests  are 
necessary  to  check  the  validity  of  this  I'^esult  for  different  flow 
geometries  and  in  different  flowfield  regions.  It  is,  of  course, 
obvious  that  if  the  model  involving  setting  the  higher-order 
correlations  to  zero  can  be  used  for  most  of  the  calculations,  it 
v/ill  lead  to  considerable  program  simplification  and  increase  in 
speed . 

The  use  of  the  ’’typical  eddy”  box  model  in  these  calculations 
has  verified  the  basic  concept  of  modeling  the  scalar  correlations 
using  a  probability  density  function  composed  of  a  set  of  Dirac 
functions  of  variable  strengths  and  position  in  the  scalar  phase 
space.  The  computations  have  also  allowed  us  to  determine  some 
shortcomings  of  the  present  model.  A  number  of  sample  modifications 
have  already  been  made  to  the  model  to  extend  its  applicability  to 
more  severe  environments ,  and  some  othei’  changes  are  currently  under 
study . 


The  current  version  of  the  [)rof’:raiTj  ir  restrict.ed  to  one-step 
forward  reactions.  In  tne  past  nionth,  we  have  developed  a  procedure 
for  handling  multi-step  reaction  systems  includinp,  the  complete 
turbulence-chemistry  interaction.  Hie  procedure  is  briefly  described 
in  the  next  section. 


III.  SF-COND-ORDER  CLOSURE  FOR  MULTI-STEP 
REACTIOIl  SYSTEMS 

The  *’t ypical  eddy box  model  has  been  most  extensively  used 
up  to  now  with  D^ree  species  (a,  3  and  y)  flov/fields.  In  principle, 
the  model  is  quite  nc^neral  and  can  be  extended  to  a  larr^  number 
of  species.  However,  the  direct  addition  of  each  new  species  in 
che  second-order  closure  system  will  add  ten  equations  for  the 
means  and  second-order  correlations  and  does  not  appear  to  be  a 
reasonable  procedure  for  complex  reaction  systems.  Actually,  such 
a  procedure  is  not  necessary. 

Consider  the  mixing  and  reaction  of  two  initially  separate 
sti^eams  a  and  S  .  These  streams  can  be  known  mixtures  of  a 
large  number  of  chemical  species.  The  streams  mix  and  react 
following  a  large  number  of  elementary  reactions  and  form  a  mixture 
of  product  species  designated  as  y  .  We  assume  that: 

1.  Ttie  composition  of  mixtures  a  and  6  is  fixed  for 
the  entire  flowfield.  There  are  no  internal 
reactions  among  the  species  in  these  two  reactant 
mixtures , 

2.  All  the  product  species  that  compose  y  are 
molecularly  mixed .  There  are  internal  reactions 
within  the  y  species  mixture.  Due  to  the  dif¬ 
ferent  rates  of  the  various  elementary  reactions 
(between  species  present  in  a  and  3  mixtures 
and  the  species  present  in  Y),the  composition  of 
the  mixture  y  varies  at  different  points  across 
the  flowfield. 

Consider  the  following  two-reaction  system  as  an  illustration 
of  the  procedure.  The  complete  set  of  reactions  that  we  will  be 
using  in  our  study  of  the  DF  chemical  laser  is  [^resented  in  Table  I. 
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L\  +  i'  -i-  I)!'  +  n 

i'\  +  n  D1-’  +  I' 

a  D,,  known 

3  iie  +  F  +  F,,  knovm  proportion;’ 
Y  Fe  +  +  F.,  +  I)  +  F  +  PP 


bi 


3 (He  +  F  +  F^) 


Y(He  +  +  F_  +  r  +  F  +  DF) 


cx(r'o)  -- 


The  interaction  between  turbulence  and  chemistry  only  has 
to  be  taken,  into  account  for  reactions  between  a  and  3  , 
a  and  y  ,  <and  6  and  y  .  Species  within  y  >'3i’e  assumed  to 
be  molocularly  mixed  anci  ,  therefoi-e,  species  correlations  do  not 
have  to  be  considered  for  the  Interiial  reactions.  For  conven¬ 
ience,  we  rewrite  the  reaction  system  as  r.ive;i  below,  taftrinp; 
each  chemical  species  with  ttie  mixtui’e  that  it  is  a  part  of: 

D,,(a)  +  F(B)  DF(y)  +  D(y) 

tl 

D.,(a)  +  F(y)  DFCy)  -t-  P(y) 

c 

D2(y)  +  F(p)  ^  PF(y)  -t-  n(Y) 


Do(y)  +  F(y) 


DI''(y)  -t  IKy)  internal  reaction 


+  P(y)  PPiY)  +  !■’(>) 

F^(y)  +  IHy)  ■’■  Pi'(Y)  +  P(y)  internal  reaction 


V(0  can  now  write  the  cVieir.i.' ;il  r-tmi’ce  tei'ma  I'o)’  tiie  component 
sreciep.  For  example, 


dP,,  (rt ) 
dt 


P.,(a)F(fi) 

P~;(a)F(Y) 

f. 


+  p:,(a)^''(p) 


+  n;(ft)F’  (y) 


dP^Cy) 

dt 


D,(y)F(B)  +  P,’{y)F'  (ItT 
i^(Y)F(Y) 


d  c  . 


Two  other  pieces  of  inforniat ion  are  r.ow  reopdred  tr^  co!:ir]ote 
this  forniulatioTi.  First,  we  will  he  selviiTf  equations  for  tfie 
second-order  eoi’re  la  t  i  ons  a  ,  and  '  Y  ^  •  t’ov;, 

D.l^(a)P^(B)  ,  ni,(a)F^(Y),  etc.,  can  be  sijnply  calculated  from 

c.  c. 

these  correlations  when  one  knows  the  cc^rripos  it  .1  or*)  of  the  nnxtures, 
and  so  the  source  terms  are  completely  determined.  Fecerri,  ;ve  have 
to  find  the  total  loss  of  5  and  P  .  In  the  above  examrle,  P 
is  beinr  lost  due  to  the  consumpt ior:i  of  both  P(P1  and  F,^(P'). 

We  determine  t)ie  term  tbxat  causes  the  greater  loss  of  P  and  use 
the  fifture.  If  at  some  point  the  consumption  of  F(P)  leads  to 
the  larf:er  loss  of  P  ,  then  one  has  to  keep  track  of  the  corres¬ 
ponding  components  F-,(P)  and  He(P)  that  now  become  part:  of  y  . 
The  composition  of  Y  has  to  be  calculated  at  eac)}  p(^]n^  across 
the  flowf iel d . 

This  procedure  operates  within  the  frarneworl-:  of  tliree  over¬ 
all  species  mixtui'es  a  ,  P  ,  and  Y  and  the  models  and  equations 
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